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Abstract 

We propose two different proofs of the fact that Oseen's vortex is the unique solution 
of the two-dimensional Navier-Stokes equation with a Dirac mass as initial vorticity. The 
first argument, due to C.E. Wayne and the second author, is based on an entropy estimate 
for the vorticity equation in self-similar variables. The second proof is new and relies on 
symmetrization techniques for parabolic equations. 

CO 

2 ' 1 Introduction 

We consider the vorticity equation associated to the two-dimensional Navier-Stokes equation, 
' namely 

d t uj{x,t) +u(x,t) ■ Vu(x,t) = Au(x,t) , x G R 2 , t>0. (1.1) 
The velocity field u(x, t) G R 2 is obtained from the vorticity m(x, t) G R via the Biot-Savart law 



u(x,t) = ±- I ( . X y ] u(y,t)dy, x G R 2 , t>0, (1.2) 
H ' 27r J R 2 \x - y\ A 



where (xi,X2)" L = (—X2,x\). It satisfies divit = and d\U2 — d2U\ = u. Equations ()1.1|) . (|1.2|) 
are invariant under the scaling transformation 

u(x,t) i ^ \ 2 uj(\x,\H) , u(x,t) i ^ Xu(Xx,X 2 t) , A > . (1.3) 

The Cauchy problem for the vorticity equation Q1.1JI is globally well-posed in the (scale 
invariant) Lebesgue space L X (R 2 ), see for instance (Sj- To include more general initial data, 
such as isolated vortices or vortex filaments, it is necessary to use larger function spaces. A 
natural candidate is the space Ai(R 2 ) of all finite real measures on R 2 , equipped with the total 
variation norm. This space contains L X (R 2 ) as a closed subspace, and its norm is invariant under 
(the spatial part of) the rescaling (jl.3|) . Moreover, A^(R 2 ) is closed with respect to the following 
weak convergence: /x n — /J, if J R2 V 9 d/x n - * J R 2 V 9 d/i for any continuous function <p : R 2 — > R 
vanishing at infinity. 
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Existence of solutions of with initial data in Ai(R?) was first proved by Cottet |1U| . 
and independently by Giga, Miyakawa and Osada ^3]) see also Kato [T5] . Uniqueness can 
be obtained by a standard Gronwall argument if the atomic part of the initial vorticity /i is 
sufficiently small but this method is bound to fail if /i contains large Dirac masses. In 

the particular case where \x = a<5o for some a E R, an explicit solution is known: 

= jG(^) , u(x,t) = -^U G (-^=) , xGR 2 , t>0, (1.4) 

where 

G(f) = i_e-l^ | „ G({) = i_|_( 1 _ e H« 1 ). {£R2 . (1 . 5) 

This self-similar solution of the two-dimensional Navier-Stokes equation is often called the Lamb- 
Oseen vortex with total circulation a. It is the unique solution with initial vorticity a5o in the 
following precise sense: 

Theorem 1.1 [H] Let T > 0, K > 0, a E R, and assume that lo G C°((0, T), L l (R 2 ) n 
L°°(R 2 )) is a solution of satisfying \\ui(-,t)\\i,i < -fT /or a// i G (0,T) and u;(-,i) — a5o as 
t^0+. Then 

w(i,t) = jG^j , xGR 2 , tG(0,T). 

Here and in the sequel, we say that u G C°((0, T), L X (R 2 ) n L°°(R 2 )) is a (mild) solution 
of (|l.lj) if the associated integral equation 

fe) = e (' 2 -' 1 ) A u;(-, ti) - r V • e(* 2 -*) A («(•, *)) dt (1.6) 

is satisfied for all < t\ < ti < T. 

The particular case covered by Theorem 11.11 is important for at least two reasons. First, it 
has deep connections with the long-time behavior of smooth solutions. Indeed, if ui(x,t) is any 
solution of (|l.lj) such that oj(-,t) G L X (R 2 ) and J* R2 w(x,t)dx = a for all t > 0, it is shown 
in ^2] that &>(•,£) converges in L X (R 2 ) to the Oseen vortex with circulation a as t — ► +cxd. 
This convergence result can in fact be deduced from Theorem II. II bv a classical renormalization 
argument, see jHj and ^1]. On the other hand, combining the standard Gronwall argument 
(which works if the initial measure has small atomic part) with Theorem 11.11 (which allows to 
treat a large Dirac mass), it is possible to prove the uniqueness of the solution of (|1.1|) for 
arbitrary data fj, G M.(R 2 ) Thus, the Cauchy problem for (jl.l|) is globally well-posed in 

the space A4(R 2 ) too. 

The proof of Theorem II .11 in |12j is based on an entropy estimate for equation (jl.ip rewritten 
in self-similar variables. For completeness, this argument will be reproduced below in its simplest 
form. The real purpose of this paper is to provide an alternative proof of Theorem ll . ll which relies 
on completely different ideas: it is based on symmetrization techniques for elliptic and parabolic 
equations which originate in the work of Talenti JH] > see [21 El • These methods ultimately rely 
on the maximum principle, whereas the original proof of ^2] uses a Lyapunov function. We find 
the comparison of both arguments very instructive, and this is why we chose to collect them in 
the present article. 

The rest of the text is organized as follows. Under the assumptions of Theorem 11.11 we 
first derive a priori estimates on the solution u>(x, t) using Gaussian bounds on the fundamental 
solution of the convection-diffusion equation, which are due to Osada in ^21- This prelimi- 
nary part is common to both approaches. Following ^21) we next introduce the self-similar 
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variables £ = x/\/t, r = log(t) and we show that the rescaled vorticity w(£,r) = e T uj(^e 2 , e T ) 
satisfies an entropy estimate which implies the conclusion of Theorem ll.il The new proof begins 
in Section |1J where we recall the definition of the symmetric nonincreasing rearrangement of a 
function. In particular, we prove that the solution of a convection-diffusion equation of the 
form (jl.ip with initial data ujq is "dominated" (in a sense to be precised) for all times by the 
solution of the heat equation with symmetrized initial data uSq . In the final section, we show 
how this key estimate provides a new and short proof of Theorem 11.11 



2 Representation and a priori estimates 

Assume that to G C°((0, T), L X (R 2 ) n L°°(R 2 )) is a solution of ifTTj) satisfying the hypotheses 
of Theorem 11.11 From we know that uj(x,t) coincides for t > with a classical solution 
of IJl.ljl in R 2 as constructed for instance in [0]. In particular uj{x,t) is smooth for t > 0. 
Since the Cauchy problem for (jl.lj) is globally well-posed in L X (R 2 ) n L°°(R 2 ) and since the L 1 
norm of uj(-,t) is nonincreasing with time, we can assume that T = +oo without loss of gen- 
erality. By assumption ||o;(-,t)||^a < K for all t > 0, hence the associated velocity field satis- 
fies t 2 \\u(-, t)\\L°° < CK for all t > 0, where C > is a universal constant, see ([S|, Theorem 2). 
The solution Lo{x,t) of (|l.lj) has an integral representation of the form 

uj(x,t) = / T u (x,t;y,s)uj(y,s)dy , x G R 2 , t>s>0, (2.1) 
Jn 2 

where T u is the fundamental solution of the convection-diffusion equation dtuj + u ■ Vuj = Aw. 
The following properties of T u are due to Osada ^7] and to Carlen and Loss jSj. 

• For any j3 £ (0, 1) there exists K\ > (depending only on K and (3) such that 

< T u (x,t;y,s) < ^j- exp(-/3 ^ ( ~_^ ) , (2.2) 

for all x,y G R 2 and t > s > 0, see (5). A similar Gaussian lower bound is also known. 

• There exists 7 G (0, 1) (depending only on K) and, for any 5 > 0, there exists K2 > 
(depending only on K and 5) such that 

|r„(x, i; y, s) - T u (x', t'; y', s')\ < K 2 (\x-xV + \t-tV /2 + \y~yV + \ssV /2 ) , (2-3) 

whenever t — s > d and t' — s' > 5, see |17j . 

• For t > s > and x, y G R 2 , 

/ r u (x,t;y,s)dx = 1 , / r u (x, t; y, s) dy = 1 . (2.4) 

J-R? JR 2 



If x, y G R 2 and i > 0, it follows from (|2.3|) that the function s 1— ► r u (ac, t; y, s) can be continu- 
ously extended up to s = 0, and that this extension (still denoted by T u ) satisfies properties (12 ,2|) 
to (|2.4|) with s = 0. Using this observation, we obtain for all x G R 2 and all t > 

w(x,t) = / T u (x,t;y,0)uj(y,s)dy 
Jr. 2 



+ 



R- 



(r u (ac, t; y, s) - r u (ac, t; y, 0)j u(y, s) dy , < s < t . 
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Since u>(-, s) is bounded in L 1 (R 2 ), it follows from ()2.3|) that the second integral in the right-hand 
side converges to zero as s goes to zero. On the other hand, since y i— > T u (x, t; y, 0) is continuous 
and vanishes at infinity, and since uj(-, s) — 1 a5o as s — > 0, we can take the limit s —* in the 
first integral and we obtain the useful representation: 

to(x,t) = aT u (x,t;0,0) , x £ R 2 , t>0. (2.5) 

This formula shows in particular that u = if a = 0. Thus, upon replacing oj(xi,X2,t) 
with — co(x2, x±, t) if needed, we can assume from now on that a > 0. To simplify the notation, 
we denote by f2(x,t) the Oseen vortex with circulation a, namely 

Q(x,t) = 4^ e_i ^ i xER 2 , t>0. (2.6) 

From (|^|) . (|2~5|) we have 

< u(x,t) < ^e^T , xGR 2 , t>0. (2.7) 

Moreover, in view of ()2.4j) . 

/ w(x,t)dx = a = / 0(x,i)dx, t>0. (2.8) 

In particular ||w(-, t)\\^i = a for all t > 0, so that we can take -fT = \a\ in the statement of 
Theorem 11.11 without loss of generality. Finally, using (|l,lj) and integrating by parts, we find 

^- [ \x\ 2 uj dx = f |x| 2 (Aa; — u ■ Vw) dx = 4a + 2 / (x ■ u)uo dx . 
d* Jn. 2 Jn? Jr 2 

Actually, the last integral vanishes identically. This can be seen by replacing u with its ex- 
pression (|1.2|) and using the symmetry properties of the Biot-Savart kernel. Thus we find 
that J r2 \x\ 2 lo dx = 4a, and in view of (|2.7j) we conclude that 

/ \x\ 2 uj(x,t)dx = Aat = \ \x\ 2 n(x,t)dx , t>0. (2.9) 
Jn 2 Jr 2 

3 Uniqueness proof using entropy estimates 

We first introduce the entropy functionals that will be used in the proof. Let / : R 2 — ► (0, +oo) 
be a C 1 function satisfying J* R2 /(£) d£ = 1. We define 



where G is as in (|1.5|) . In kinetic theory, is the relative Boltzmann entropy of the distribu- 

tion function / with respect to the Gaussian G. In information theory, H(f) is often called the 
relative Kullback entropy of / with respect to G, and 1(f) the relative Fisher information of / 
with respect to G. The entropy H(f) satisfies the following bounds (see for instance 0]) 

\\\f-G\\ 2 Ll < H(f) < 1(f) , (3.2) 

where the first inequality is the Csiszar-Kullback inequality and the second bound is the Stam- 
Gross logarithmic Sobolev inequality. It follows in particular from (|3.2j) that H(f) > 0, 
and H(f) = if and only if / = G. 
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Now, let u>(x, t) be a solution of satisfying the hypotheses of Theorem ll.il and let u(x, t) 
be the associated velocity field. As in the previous section, we assume without loss of gener- 
ality that T = +00 and that a = J R2 Ld(x,t) dx > 0. Following ^21) we define the rescaled 
vorticity u)(£,r) and the rescaled velocity field ?;(£, r) by 

w(£,t) = e T uj(£ei,e T ) , v(£,r) = eiu(£ei,e r ) , £ G R 2 , rGR. (3.3) 

The equation satisfied by u>(£,r) reads 

<9 r w + (v ■ Vz)w = A^w + -(£ • V|)w + u; . (3.4) 

Since divf = and d\V2 — 02^1 = w, the rescaled velocity v is obtained from the rescaled 
vorticity w via the usual Biot-Savart law. 

Our aim is to show that u>(£,r) = aG(^) for all £ G R 2 and all r G R. To do that, we 
consider the relative entropy h(r) = H(w(-,T)/a), where H is defined in ()3.1j) . Since by (12.71) 
we have < w(£,t) < Kiae'^ 2 ^, there exists K 3 > such that < h{r) < K 3 for all r G R. 
On the other hand, using (J3.4|) and integrating by parts, it is not difficult to verify that 

^H(^l) = -l(^l) , rGR, 
dr V a J \ a J 

see [TSJ. In view of the second inequality in ()3.2|) . this implies that h'{r) < —h(r) for all r G R. 
Thus 

/i(t) < e" {r " To) /i(r ) < ^ 3 e" {r - To) whenever r > r . 

Letting To — > —00, we obtain /i(t) = for all r G R, hence «;(£, r) = aG(£) for all £ G R 2 and 
all r G R. This is the desired result. □ 



4 Symmetric rearrangements and applications 

In this section, we recall the definition of the symmetric nonincreasing rearrangement of a 
function / : R^ — ► R, and we establish two results that will be used in the final section. We 
refer to pQ, 0, |16| for more information on rearrangements and their applications. 

Let / : R^ — ► R be a measurable function. We assume that / vanishes at infinity in 
the sense that meas({x G H N \ \ f(x)\ > £}) < 00 for all t > 0. We define the distribution 
function fif : [0, +00) — ► [0, +00] by 

H f (t) = meas({x G R^ | \f(x)\ > t}) , t > , 

and the nonincreasing rearrangement /* : [0, +00) — * [0, +00] by 

f*(s) = sup{i > I n f (t) > s} , s > . 

The symmetric nonincreasing rearrangement : R^ — > [0, +00] is then defined by 

f*(x) = f*(c N \x\ N ) , x G H N , 

where c N = ir N / 2 /T(f + 1) is the measure of the unit ball in R^. By construction, is 
radially symmetric and nonincreasing (along rays). Moreover, is lower semicontinuous (hence 
measurable) and meas({x G H N | /*(x) > t}) = meas({x G R^ | \f(x)\ > t}) for all t > 0. As a 
consequence, if / G LP(R N ) for some p G [1, 00], then f* G ^(R^) and = \\ fill*. More 

generally, one has ||/* — < ||/ — g\\jj> for all /, g G L P (R N ). Note also that /* = \f\^, 

and that is a continuous function if / is continuous. 

We next introduce a partial order on integrable functions which is based on rearrangements. 
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Definition 4.1 Let f,g : H N — ► R 6e integrable functions. We say that f is dominated by g if 

[ f*{x)dx < I g#(x)dx forallR>0, (4.1) 
Jb r J b r 

where Br = {i£ H n \\x\ < R}. In this case, we write f ■< g. 
It can be shown that / ^ g if and only if 



<S>(\f(x)\)dx < / $(\g(x)\)dx 
r n Jr n 

for all convex functions : [0, +00) — > [0, +oo) with $(0) = 0, see pQ. In particular, if / ^ g 
and g G L 1 (R Ar ) n ^(R^), then ||/|| LP < || 5 || LP for any p > 1. 

The uniqueness proof in Section El will be based on two properties of the domination relation 
which we now describe. These results are rather standard, and the proofs will be outlined for 
completeness only. 

Proposition 4.2 Let f,g : H N — ► [0, +oo) be continuous and integrable functions satisfying: 

a) f < g; 

b) g = g*; 

c ) Jrn f( x ) dx = Irn 9{x) dx; 

d) / R]V \x\ N f{x)dx = f RN \x\ g{x) dx < co. 
Then f = g. 

Remark 4.3 As is clear from the proof, the same result holds if \x\ N is replaced by in 
assumption d), where <p is any continuous, positive, and strictly increasing function on [0, +oo). 

Proof: We first observe that 

x\ N f*(x)dx< [ \x\ N f(x)dx, (4.2) 

R N JR N 

and that equality holds in Q4.2JI if and only if / = f*. Indeed, this is obvious if / is the 
characteristic function of a bounded open set A C R^, because then f & = 1b r where cjyR N = 
meas(A). The general case follows using the "layer cake representation" |16j 



f(x) = J l {f>t} (x)dt, f*{x) = J lf f>t} (x)dt, 



together with Fubini's theorem. 

Now, let h = g* — f*, where f*,g* are the nonincreasing rearrangements of /, g. For all R > 0, 
we have 

p rc N R N r- i rc N R N 

/ f#{x)dx = / f*(s)ds, / \x\ N f*{x)dx = — sf*{s)ds . 

J b r Jo J b r °n Jo 



Using the first equality, we can rewrite assumptions a) and c) as follows: 
a') J r h(s) ds>0 for all r > 0, 
c') / °° h(s) ds = 0. 

In view of (|4.2|) and the second equality, assumptions b) and d) also imply 
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d') Jo°° s ^( s ) ds > 0, with equality if and only if / = 



Let W(r) = Jq h(s) ds = - J r °° ds > 0. Clearly, rW(r) < J r °° s|/i(s)| ds -► as r -c +oo, 
hence 

r=oo 

= . 



sh(s)ds + H(s)ds = / (s/i(s) + H(s)) ds = rH{r) 



r=0 



By a') and d'), this implies that J °° sh(s) ds = (hence / = /*) and 7i = (hence we 
have = 5* = g). This concludes the proof. □ 

Assume now that / : x [0, +oo) -> R is a solution of the convection-diffusion equation 

dtf(x, t) + U(x, t) ■ V/(x, i) = A/(x, t) , xeR N , t>0, (4.3) 

where U : K N x [0,+oo) is a smooth, divergence-free vector field which is bounded 

together with all its derivatives. Assume moreover that the initial condition fo(x) = f(x,0) is 
continuous and decays rapidly at infinity. 

Proposition 4.4 Under the assumptions above, the solution of &4-3Q with initial data /o satisfies 

f(;t) 1 e tA f* , forallt>0. 
Proof: The proof uses two main ingredients: 

i) If / 9 and g = g#, then e tA f < e tA g for all t > 0, see ([3], Proposition 3). This is because 
the integral kernel of the heat semigroup is positive, radially symmetric, and decreasing along 
rays (in other words, it coincides with its own rearrangement). 

ii) Let S(t\,t2) be the evolution operator associated to the linear nonautonomous transport 
equation 

d t f{x,t) + U(x,t) -Vf(x,t) = , xe~R N , t>0, 

so that f(-,ti) = S(ti,t2)f(-,t2). Since U(x,t) is globally Lipschitz and divergence- free, one 
has S(ti,t2)f = f o ^(ti,t2)~ l , where ^(ii^) : R- N R N is a measure preserving transfor- 
mation. It follows immediately that (»S(ii, £2)/)* = 

To prove the claim, we use an approximation procedure in the spirit of Trotter's formula, 
see (|19|. Section 11. A). Fix t > 0. Given n G N*, we define finite sequences {f£}k=o,...,n 
and {g%}k=o,..., n inductively by 

f n _ q ( ( k + l ) t kt\ ± A 

Jk+i ~ H n '^J e ' Jo -Jo, 

„n „— A„n „n f# 

g k+ i - en g k , g - Jo ■ 

Using properties i) and ii) above, it is easy to verify by induction over k that fj} X g^ for 
all k G {0, ...,n}. In particular, ffi ^ g™. But g™ = e* A /*, and a classical result shows 
that /" — ► f(-,t) in L 1 (R N ) as n — > oo, so that (/")* — > /*(•, i) also. Thus passing to the limit 
we obtain /(-,t) ^ e* A / # . □ 



5 Uniqueness proof using symmetrization 

Let again w(ar, i) be a solution of (|1.1|) satisfying the hypotheses of Theorem ll.H and let u(x, t) be 
the associated velocity field. As usual we assume that T = +oo and that a = J R2 u(x, t) dx > 0. 
For any t > 0, let u>&(x, t) be the symmetric nonincreasing rearrangement of u(x, t) with respect 
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to the spatial variable x E R 2 . Using (|2.7|) and the fact that the rearrangement is order 
preserving, we find 

< LO*(x,t) < e _/3 ^ , xGR 2 , t>0. 

Since J R2 u>#(x, t) dx = J R2 u>(x, t) dx = a, we deduce that u>#(-,t) — a^o as i — > 0+. 

Fix t > s > 0. Applying Proposition 14.41 with N = 2, f(x,t') = uj(x,t' + s), and U(x,t') = 
u(x, t! + s), we obtain 

Lo(-,t) < e^ A io*(-,s) . 

As s — > 0+, the right-hand side converges pointwise (hence also in L X (R 2 ) by the dominated 
convergence theorem) to the Oseen vortex £l(-,t) defined in (|2,6|) . It follows that 

w(-,t) r< n(-,t) , for alH > . (5.1) 

We now fix t > and apply Proposition 14.21 with N = 2, / = u>(-,t), and 5 = f2(-,i). All 
assumptions are satisfied due to ()2.8I) . (|2.9|) . and ()5.1() . We conclude that u>(x,t) = £l(x,t) for 
all x S R 2 , which is the desired result. □ 
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